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Based on previous observations of glaze ice accretion on aircraft surfaces, a multlzone mode! with distinct 
zones of different surface roughness is demonstrated. The use of surface roughness In the LEWICE Ice accretion 
prediction code Is examined. It was found that roughness Is used in two ways: I) to determine the laminar to 
turbulent boundary-layer transition location and 2) to calculate the convective turbulent heat-transfer coefficient. 
A two-zone version of the multizone model is implemented in the LEWICE code, and compared with experimental 
convective heat-transfer coefficient and ice accretion results. The analysis of the boundary-layer transition, 
surface roughness, and viscous flowfleld effects significantly increased the accuracy in predicting heat-transfer 
coefficients. The multizone model was found to significantly improve the ice accretion prediction for the cases 
compared. 


I. Introduction 

A. Background 

I n rime ice accretion, which occurs at cold temperatures, 
there is adequate convective heat transfer to rapidly freeze 
the impinging supercooled water droplets. The rapid freezing 
results in an opaque white ice due to entrapped air bubbles. 
Because of the rapid freezing, rime ice shapes are determined 
solely by the droplet impingement behavior; efforts to model 
rime ice accretion have been relatively successful. 

In contrast, current efforts to analytically model glaze ice 
accretion are hampered by insufficient knowledge of the ac- 
creting ice surface roughness. In glaze icing, which normally 
occurs at temperatures near freezing or at high liquid water 
contents, there is insufficient convective heat transfer to re- 
move all of the latent heat of freezing from the impinging 
supercooled water droplets. The water freezes slowly, re- 
sulting in a strong clear ice structure. The local ice accretion 
rate is controlled by the ability to remove latent heat and thus 
the local convective heat transfer. The droplet impingement 
behavior also plays an important role in determining where 
on the surface there is sufficient impingement that the surface 
is wet. 

The local convective heat transfer from a surface is known 
to be strongly dependent on the ice surface roughness. 1 Be- 
cause of the importance of heat transfer on the ice accretion 
rate, surface roughness becomes an important factor in mod- 
eling glaze ice accretion. Current analytical models such as 
LEWICE generally assume that the surface roughness is uni- 
form and the effective sand grain roughness k, is used as an 
input parameter in the code. 2 - 3 The magnitude of k t> the 
roughness parameter, is normally determined empirically by 
comparison of predicted and experimental ice accretions. The 
erratic performance of glaze ice accretion models and the 
empirical manner in which the surface roughness is treated 
indicate the need for a more deterministic treatment of the 
surface roughness. 

B. Multizone Model 

In prior experiments, detailed observations of accreting ice 
surface roughness were made at several icing facilities/ 6 As 
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a result of the detailed photographic analysis of accreting ice 
surfaces in these facilities, distinct regions were identified each 
having a characteristic roughness and identifiable boundaries. 
Based on these observations, a multizone mode], in which the 
accreting ice surface is divided into two or more discrete zones 
with varying surface roughness and water behavior, was pro- 
posed. 

/. Surface Roughness Zones 

A total of three roughness zones were identified. A typical 
ice accretion showing the three zones is shown in Fig. 1. 

1) Smooth zone'. Close to the stagnation point, the surface 
was observed to be uniformly wet with a thin film of water 
at warm temperatures. In this “smooth zone,” the surface 
was smooth with no distinctively visible roughness. 

2) Rough zone\ At some point downstream, there was a 
sudden transition to a significantly rougher zone, where there 
appeared to be insufficient water to maintain a uniform film, 
and the water tended to coalesce into the water beads first 
observed by Olsen and Walker. 78 In this “rough zone,” the 
ice accretion rate was observed to be enhanced compared to 
the smooth zone. This is thought to be due to increased heat 
transfer resulting from the greater surface roughness. 

In certain glaze ice accretions, particularly at high liquid 
water contents, some of the roughness elements within all or 
part of the rough zone were observed to grow into distinct 
protuberances. With higher collection efficiency and en- 
hanced heat transfer, these protuberances tended to grow 
rapidly resulting in horned ice accretions. 

3) Runback or rime feather zone: Behind the impingement 
region, aft of the rough zone, ice was observed to accrete by 
two different processes, depending on temperature. At warm 
temperatures, a runback zone was observed aft of the rough 
zone, which was characterized by areas of ice interspersed 
with un-iced surface. The surface water was observed to ini- 
tially runback and then stagnate at the point of flow sepa- 
ration, where this water slowly froze as rivulets or as large 
coalesced water cells. At cold temperatures, rime feathers 
were sometimes observed to grow in the region aft of the 
primary accretion. The feathers were dry ice accretions, which 
propagated in the local upstream direction. 

2. Roughness Transition Location 

The angular position of the boundary between the smooth 
and rough zones has been experimentally observed as a func- 
tion of time for cylinders. A typical example is shown in Fig. 

2 for a 1 in. diameter cylinder. This data was obtained using 
the technique discussed in Sec. IVA. The initial location of 
roughness transition is believed to be caused by, and to co- 
incide with, the location of the boundary-layer transition from 
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Fig. 1 Typical glaze ice shape for a 2-min exposure on a 1-in. cylinder 
showing the three distinct roughness zones ( T = — 4.5°C, V = 150 
kt, liquid water content = 1.0 g/m\ and median volumetric diameter 
= 30 pm). 
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Fig. 2 The angular position vs time of smooth-rough transition lo- 
cation (T = -7*C, V = 125 kt, liquid water content = 0.8 g/m 5 , and 
median volumetric diameter = 12 pm). 


laminar to turbulent. Experimentally, the rough/smooth tran- 
sition was observed to propagate with time towards the stag- 
nation region for a wide variety of icing conditions. The effect 
of icing cloud parameters on the transition location has been 
studied in an attempt to identify the underlying physical mech- 
anisms that cause the rough surface to develop. It was found 
that higher velocity tends to move the initial location of the 
transition closer to the stagnation point. Warmer tempera- 
tures and higher liquid water contents tend to move the tran- 
sition point more quickly towards the stagnation region. 

Observations that the smooth/rough transition depends on 
Reynolds number indicate that the initial transition in surface 
roughness is caused by boundary-layer transition. In these 
cases, the laminar/turbulent boundary-layer transition point 
will also be the initial smooth/rough transition point. In the 
laminar region, the heat transfer is low enough that all of the 
latent heat cannot be removed and the surface is coated by 
a uniform water film. However, the enhanced heat transfer 
in the turbulent region can cause sufficient freezing to partially 
dry the surface and cause bead formation, which results in 
the increased roughness and the increased heat transfer ob- 
served in the rough region. 

The influence of surface water flow on the smooth/rough 
transition location also indicates that dynamic effects are im- 
portant. It is thought that the dynamic effects are caused by 
bead formation at the interface between the smooth and rough 
surface zones as shown in Fig. 3. The formation of these beads 
at the interface causes enhanced heat transfer within the rough 
zone, which tends to freeze out the downstream beads and 
dry the surface. As the surface dries, beads begin to form 
further upstream and the transition point will propagate to- 
wards the stagnation region as observed experimentally. As 
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Fig. 3 Schematic representation of bead formation at the smooth-to- 
rough transition. 


roughness moves upstream, the boundary layer is tripped and 
will move with the roughness transition. By increasing the 
surface water flux, the rate of formation and growth of the 
interfacial beads is increased. This causes the observed faster 
upstream propagation of the transition point with increased 
surface water flux. 

Based on the experimental observation of the smooth/rough 
transition location, it was determined that for some cases, the 
initial smooth/rough transition occurs on the laminar/turbu- 
lent boundary-layer transition point. Therefore, for the pur- 
poses of ice accretion prediction, it is hypothesized that the 
smooth/rough transition location can be assumed to coincide 
with the laminar/turbulent boundary-layer transition. 

3. Model Description 

In contrast to the conventional techniques that assume uni- 
form roughness over the entire ice accretion, the multizone 
model is divided into two or more discrete zones with varying 
surface roughness and surface water behavior to be consistent 
with experimental observations. 

In the simplest version of the model, the surface is divided 
into two zones, the “smooth” zone and the “rough” zone. In 
the “smooth” zone, corresponding to the smooth region cen- 
tered about the stagnation line, the surface is uniformly wet, 
with thin water film runback. The heat transfer is that for a 
laminar boundary layer and a Messinger 9 type water runback 
model used in the original LEWICE code appears to be valid. 

In the “rough” zone, surface tension effects are important 
and the characteristic water beads or roughness elements ap- 
pear. Here, the heat transfer is enhanced due to increased 
roughness and the experimental results have indicated no water 
runback in this region for some cases. 5-8 

II. Implementing the Multizone Model Through 
Boundary-Layer Transition 

Based on the hypothesis that the smooth/rough transition 
location coincides with the laminar/turbulent boundary-layer 
transition, a simple more physically realistic multizone model 
was implemented in the LEWICE code. Basically, the smooth 
and rough zones were considered to coincide with the regions 
of laminar and turbulent boundary layer, which are treated 
separately when calculating heat transfer in LEWICE. In or- 
der to implement the multizone model effectively, the use of 
surface roughness in the LEWICE code was examined. The 
multizone model was then implemented in the LEWICE code 
through boundary-layer transition in order to evaluate the 
effectiveness of the model. 

A. Use of Roughness In LEWICE 

In the original LEWICE code, the only roughness param- 
eter is the equivalent sand grain roughness height k s . This 
parameter, which is one of the most important factors in 
determining the ice shape, 7 8 is used in two ways. One use is 
to determine the location of boundary-layer transition from 
laminar to turbulent, and the other is to calculate the mag- 
nitude of heat transfer in the turbulent region. 

The boundary-layer transition is determined in the laminar 
region with reference to a critical Reynolds number based on 
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Fig. 4 Effect of k, in the laminar region on the transition location. 
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Fig. 5 Effect of k. In the turbulent region on the magnitude of the 
heat-transfer coefficient (Note: Normalized HTC - HTC/VNu). 


the roughness element height and local velocity. This Rey- 
nolds number is 



where k , is the roughness element height in the laminar region, 
v k is the air velocity at distance k , from the surface, and v is 
the kinematic viscosity of the air. When this local Reynolds 
number exceeds a critical value of 600, the boundary layer is 
transitioned to turbulent, 10 and any region downstream of this 
point is considered to have turbulent boundary layer. The 
effect of varying k t in this region on the transition location is 
demonstrated in Fig. 4 by running LEWICE on a 0.15-m- 
diam cylinder at a velocity of 25 m/s and temperature of 0°C. 
It is noted that no transition occurs until a value of k, above 
0.5 mm. At higher values, the transition location propagated 
towards the stagnation region with increasing k s . 

In the laminar region, the heat transfer is independent of 
the roughness element height, assuming the roughness height 
is less than the boundary-layer thickness. In regions with a 
turbulent boundary layer, the roughness element height is 
used in another way. Here, roughness plays an important role 
in the heat transfer. Calculation of the heat-transfer coeffi- 
cient in the turbulent region involves a complicated set of 
equations based on flow parameters, roughness, and the ge- 
ometry of the accreting body, including the roughness element 
height. 3 Because of the complexity of these equations, it is 
very difficult to determine the effects of the equivalent sand 
grain roughness k t on the turbulent heat transfer analytically. 
However, these effects of k, on heat transfer were observed 
numerically using the LEWICE code. An example is shown 
in Fig. 5 where the convective heat-transfer coefficient (HTC) 
normalized by the square root of the Nussult number (i.e., 
normalized HTC = HTC/VNu) is plotted against /c Jf for a 
0.15-m-diam cylinder at a velocity of 33 m/s and temperature 


of 0°C. It was found that increasing k s significantly increased 
the heat transfer in the turbulent region. 

Because of the manner in which roughness is treated in the 
LEWICE code, it is possible to separate the surface roughness 
in the laminar and turbulent regions. The laminar roughness 
value is then used to determine the laminar to turbulent 
boundary-layer transition location. The turbulent roughness 
is used to calculate the heat-transfer coefficient in the tur- 
bulent region. 

B. Implementation of the Multizone Model 

In order to evaluate the effectiveness of the multizone model, 
a simple two-zone version was implemented in the LEWICE 
code. For the first time step, different roughness are used for 
the laminar and the turbulent regions. In subsequent time 
steps, the boundary-layer transition point is made to propa- 
gate towards the stagnation line in a manner consistent with 
experimental observations of the smooth/rough transition. In 
order to avoid errors due to improper boundary-layer tran- 
sition placement in these initial evaluations, the transition 
location was input to the code from the experimental mea- 
surements of smooth/rough transition migration such as those 
shown in Fig. 2. It should be noted that a deterministic pre- 
diction of the smooth/rough transition migration valid for all 
geometries and based on nominal code inputs is necessary for 
a full implementation of the multizone model. 

For the initial time step, the user can input two different 
roughness element sizes, one for the laminar region k t and 
another for the turbulent region k r As explained in Sec. IIB, 
the roughness element height is used in LEWICE for different 
purposes in the two regions. By separating the two, it is pos- 
sible to control the transition location and the turbulent heat 
transfer independently. This method more closely emulates 
the physical situation where two different roughness element 
heights have been observed in the smooth and rough regions. 
In the laminar region, a roughness height k, corresponds to 
the roughness of the uniform water film found in the exper- 
imental observations. In the turbulent region, the roughness 
height k, corresponds to the roughness size observed in the 
rough zone. This method is valid only for the first time step 
where the transition location is determined mainly from the 
boundary-layer instability and the dynamic effects of surface 
water are negligible. 

For subsequent time steps, the angular locations of the 
laminar to turbulent transition point as a function of time are 
specified based on the experimental data. The roughness ele- 
ment height k r for the turbulent region still needs to be spec- 
ified. For some cases, however, the experimental results have 
indicated that freezing fraction is unity (i.e., all impinging 
water freezes) in this region. For these cases, the ice shape 
is not dependent on the roughness element height that user 
chooses, since it, is large enough to ensure that there is enough 
heat transfer to remove all of the latent heat. 


III. Ice Shape Comparisons on Cylinders 

A. Comparison with Original LEWICE 

In order to investigate the discrepancy that had been ob- 
served between the experimentally observed roughness ele- 
ment sizes and the recommended roughness element size in 
LEWICE, 4 several original LEWICE runs were compared to 
the experimental heat-transfer coefficient data of Achenbach 
for rough 0.15-m-diam cylinders. 1 The measurements provide 
convective heat transfer coefficient data for cylinders with 
known surface roughness sizes at various Reynolds numbers. 
In the following figures, Vindicates the angular position along 
the cylinder, where 0 deg is the stagnation point. The vertical 
axis indicates the heat transfer coefficients normalized by the 
square root of the Nussult number. 

A typical heat transfer comparison is shown in Fig. 6 for a 
case with Re = 1.27 x 10 5 and a moderate roughness element 
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Fig. 6 LEWICE predicted and Achenbach measured heat-transfer 
coefficients for a cylinder with equivalent sand grain roughness k , = 
0.115 mm and Re = 5.9 x 10* (Note: Normalized HTC - 

HTCIVNu). 



8 (d*g) 

Fig. 7 LEWICE predicted and Achenbach measured heat-transfer 
coefficients for a cylinder with equivalent sand grain roughness k t ~ 
0.45 mm and Re - 1.27 x 10 5 (Note: Normalized HTC ~ 

HTCtVNu ). 


size k t of 0.45 mm. In this case, the heat transfer is well- 
predicted in the laminar region. However, where the transi- 
tion to a turbulent boundary layer is indicated by a rise in the 
Achenbach heat transfer data at 75 deg, LEWICE predicts 
transition further forward at 57 deg. This discrepancy can be 
explained by one of the following hypotheses. One is that the 
roughness element height k, used to calculate the critical tran- 
sition Reynolds number is not correct. The other is that the 
transition model, based on the critical roughness element Rey- 
nolds number of 600, is not valid for this case. 

The comparisons shown in Fig. 6 also indicate a difference 
in the magnitude of heat transfer predicted in the turbulent 
region. LEWICE can be seen to clearly overpredict the heat 
transfer in this region. This result indicates the need to further 
investigate the effect of roughness on heat transfer in the 
turbulent region. 

Another comparison is shown in Fig. 7 for a case with Re 
= 5.9 x 10 3 and a small roughness element height k, of 0.1 15 
mm. Here, LEWICE does not predict transition at all, where 
the experimental result indicates a transition at 63 deg. For 
this case of very small roughness element size, the transition 
model is clearly not valid, since it fails to predict transition 
at all, indicating a limit of the critical transition Reynolds 
number theory. 

B. Effect of External Flowfleld Model 

Because the local convective heat transfer coefficient is 
strongly dependent on the external flow velocity on the sur- 
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Fig. 8 A comparison of normalized heat-transfer coefficients calcu- 
lated viscous and inviscid flowfields (Re - 7.3 x 10 4 and = 1.35 
mm). 


face, errors in the velocity field will manifest as errors in the 
convective heat transfer coefficient. This can be seen in Fig. 
6. Because LEWICE uses an inviscid flow model, it cannot 
accurately calculate the flowfield in the separation region be- 
hind the cylinder. This is thought to be the cause of the large 
discrepancy between LEWICE and the experimental results 
in Fig. 6 for angles greater than 135 deg. 

The inability of LEWICE to accurately calculate heat trans- 
fer in separation and recirculation zones may be significant 
for glaze ice modeling where such regions are common behind 
“horn-shaped'’ ice accretions. It should be noted, however, 
that there is very little droplet impingement in the recircu- 
lation regions and that the heat transfer discrepancies will be 
most apparent in ice accretions where there is surface water 
runback into the recirculation region. 

In order to verify that it is possible to predict heat transfer 
more accurately using the LEWICE methodology with a more 
accurate velocity field, the potential flowfield was replaced 
with a well-documented empirical formula for cylinders. 11 The 
heat transfer was calculated using the empirical flowfield and 
an example is shown in Fig. 8. In this comparison, the critical 
Reynolds number was set at 975 to match the boundary-layer 
transition location based on the known roughness element 
height of 1.35 mm. It can be seen that the heat-transfer pre- 
diction was significantly improved by the use of the viscous 
flowfield and a more accurate prediction of boundary-layer 
transition location. 


IV. Comparison with Experimental Ice 
Accretion Data 

A. Experimental Setup 

In order to obtain roughness transition locations as a func- 
tion of time on well-defined ice shapes for comparison with 
the multizone model, detailed observations of accreting ice 
shapes on cylinders at various cloud conditions were made at 
the B. F. Goodrich Ice Protection Research Tunnel. A series 
of l-in.-diam cylinders were observed at a freestream velocity 
of 125 kt. The liquid water content (LWC) ranged from 0.4- 
1.5 g/m 3 with median volumetric diameter (MVD) of 11-40 
/xm. The tunnel total air temperatures varied from 0°C to 
- 15°C. 

The photographic setup used in the tunnel is shown in Fig. 
9. The cylinders horizontally spanned the test section. The 
test section walls had heated glass windows to provide pho- 
tographic access. Two charge-coupled device (CCD) high- 
magnification microvideo cameras were used to obtain a graz- 
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Fig. 9 Schematic diagram of the photographic setup in the icing wind 
tunnel. 
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Fig. 10 Schematic diagram of the fan laser beam setup In the Icing 
wind tunnel. 


ing angle view of the ice accretion. The cameras were focused 
at the central region of the tunnel where the cloud was most 
uniform. 

A laser light sheet was used to illuminate a plane perpen- 
dicular to the test cylinders to accurately identify the plane 
of focus (Fig. 10). With the laser sheet, it was possible to 
record on video and still cameras the details of accreting ice 
shapes and roughness at a single spanwise location. Because 
the rougher surface in the rough region increased the internal 
reflection, the ice shape appeared brighter in the rough re- 
gion, which allowed an accurate identification of smooth to 
rough transition point. The external light inside the tunnel 
was turned off at a specific time interval, usually 20 s, pro- 
viding enhanced details in the laser illuminated ice shape. 

B. Comparison of Actual and Predicted Ice Shapes 

An example of actual and predicted ice shapes using the 
multizone LEWICE code and the original LEWICE code 
values for a 4-min ice accretion at 20-s intervals is shown in 
Fig. 11. The icing cloud condition was -7°C, 125 kt, LWC 
of 0.8 g/m 3 , and a MVD of 12 /im. The angular position of 
transition was measured from the video recordings at 20-s 
intervals and is shown in Fig. 2. Characteristically, the tran- 
sition was observed to migrate towards the stagnation region 
with time. For the original LEWICE case which used a 
recommended 3 roughness height value of 0.04 mm, a rela- 
tively thin conformal ice accretion was predicted. However, 
for the multizone case, the smooth and rough ice accretion 
zones significantly improved the prediction of the overall ice 
shape, showing a good correlation between the experimental 
ice accretion and predicted ice accretion using the multizone 
version of LEWICE. Although additional validation cases 
must be run, the initial improvement of the multizone pre- 
diction is encouraging. 



a) Original LEWICE prediction 




c) Multlzone LEWICE prediction 

Fig. II Comparison of original and multlzone LEWICE predictions 
with the experimental ice shapes for a I -in. cylinder at 20-s intervals 
(T - -7*, V = 125 kt, liquid water content = 0.8 g/m 3 , and median 
volumetric diameter = 12 /un). 


V. Conclusions 

The investigation of the heat transfer on accreting ice sur- 
faces has resulted in the following conclusions; 

1) The use of roughness element height in LEWICE has 
been investigated. The roughness influences the location of 
the laminar to turbulent boundary- layer transition and the 
convective heat transfer in the turbulent bound ary- layer re- 
gion. The convective heat transfer in the laminar region is 
independent of the roughness element height. 

2) Increasing the roughness height in the laminar region 
was observed to move the laminar to turbulent boundary- 
layer transition location closer to the stagnation region. In- 
creasing the roughness height in the turbulent region in- 
creased the turbulent heat transfer. 

3) The heat transfer predicted by the multi-zone version 
of LEWICE was compared with known data on cylinders. It 
was found that by using a viscous fiowfield and appropriate 
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boundary-layer transition criteria, it was possible to predict 
the heat transfer quite accurately. This indicated the impor- 
tance of having a viscous flowfield model in regions of sep- 
aration and recirculation. 

4) A laser fan beam technique was used to obtain high- 
resolution profiles of ice accretions and to highlight the smooth- 
to-rough surface transition location. This technique signifi- 
cantly improved the ease and accuracy of ice accretion pho- 
tography. 

5) An initial multi-zone version of the multizone model 
was implemented in the LEWICE ice accretion prediction 
code. For the first time step, different roughness heights for 
laminar and turbulent regions were used. For the subsequent 
time steps, an experimentally observed smooth-to-rough tran- 
sition location was used to impose laminar to turbulent bound- 
ary-layer transition. The experimental observations were used 
to verify the concept of the multizone model. The multizone 
model significantly improved the prediction of the glaze ice 
accretion. However, further work is required to validate the 
approach and to develop a deterministic prediction of the 
rough- to-smooth transition dynamics. 
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